Collective states in highly symmetric atomic configurations, and single-photon traps 



Hanno Hammer* 

Department of Mathematics 
University of Manchester Institute of Science and Technology ( UMIST) 
P.O. Box 88 
Manchester M60 1QD 
United Kingdom 
(Dated: February 1, 2008) 

We study correlated states in a circular and linear-chain configuration of identical two-level atoms containing 
the energy of a single quasi-resonant photon in the form of a collective excitation, where the collective behaviour 
is mediated by exchange of transverse photons between the atoms. For a circular atomic configuration contain- 
ing N atoms, the collective energy eigenstates can be determined by group-theoretical means, making use of the 
fact that the configuration possesses a cyclic symmetry group Zjv- For these circular configurations the carrier 
spaces of the various irreducible representations of the symmetry group are at most two-dimensional, so that the 
effective Hamiltonian on the radiationless subspace of the system can be diagonalized analytically. As a con- 
sequence, the radiationless energy eigenstates carry a Zjv quantum number p = 0, 1, ... ,N which is analogous 
to the angular momentum quantum number I = 0, 1, . . ., carried by particles propagating in a central potential, 
such as a hydrogen-like system. Just as the hydrogen ^-states are the only electronic wave functions which can 
occupy the central region of the Coulomb potential, the quasi-particle corresponding to a collective excitation 
of the circular atomic sample can occupy the central atom only for vanishing Z;v quantum number p. When a 
central atom is present, the p = state splits into two, showing level-crossing at certain radii; in the regions be- 
tween these radii, damped oscillations between two "extreme" p = states occur, where the excitation occupies 
either the outer atoms or the central atom only. For large numbers of atoms in a maximally subradiant state, a 
critical interatomic distance of X/2 emerges both in the linear-chain and the circular configuration of atoms. The 
spontaneous decay rate of the linear configuration exhibits a jump-like "critical" behaviour for next-neighbour 
distances close to a half-wavelength. Furthermore, both the linear-chain and the circular configuration exhibit 
exponential photon trapping once the next-neighbour distance becomes less than a half-wavelength, with the 
suppression of spontaneous decay being particularly pronounced in the circular system. In this way, circular 
configurations containing sufficiently many atoms may be natural candidates for single-photon traps. 
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I. INTRODUCTION 



In the classical paper by Dicke [1], super- and subradiance in a collection of two-level atoms was studied in the near- field 
regime from a theoretical point of view. Within the same near-field limit, the topic has subsequently been reviewed by a couple 
of authors [2, 3], with a comprehensive study of subradiance being given in the series of papers [4-7]. While Wigner functions, 
squeezing properties and decoherence of collective states in the near- field regime have been presented in [3, 8, 9], triggering 
of sub- and superradiant states was shown to be possible in [10]. An experimental observation of super- and subradiance was 
reported in [1 1, 12]. The near- field limit utilized in these examples is called the "small-sample approximation", the term deriving 
from the fact that the atoms are assumed to be so close together that they are all subject to the same phase of the radiation field. 
In this approximation, the interaction Hamiltonian is independent of the spatial location of the atomic constituents, similar to a 
long-wavelength approximation. 

In contrast, our work presented here focuses on super-and subradiance in highly symmetric atomic systems with arbitrarily 
large interatomic distances. While in the first part of this report we deal with the theory of simply-excited correlated states in 
an arbitrary sample of atoms, in the second part we shall mainly be concerned with subradiance and, in particular, with the 
radiation-trapping capability, of circular and linear-chain atomic configurations. This radiation trapping is a consequence of 
collective states of the atomic sample which can have extremely small spontaneous decay rates as compared to the single-atomic 
decay rate; these states therefore are natural candidates for single-photon traps. 
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As mentioned above we are mainly interested in the decay of simply-excited states of the atomic sample; these states may 
be conceptualized as follows: Consider N identical two-level atoms with infinite mutual distances such that precisely one of 
the atoms is excited, all others being in the ground state. If we now think of adiabatically decreasing the distances, the atoms 
will begin to interact with each other via their coupling to the radiation field, and hence the excitation, which previously was 
localized at one atom only, will distribute over the whole sample. The resulting states therefore will be superpositions of the 
excited levels in the atoms in such a way that the atomic sample still contains the energy of one (more or less resonant) photon, 
but this energy is now delocalized over the whole sample. Moreover, as a consequence of this derealization, the atoms will be 
strongly entangled with each other as well as with the radiation field. Some of these correlated states turn out to be capable of 
very strong suppression of spontaneous decay; this radiation trapping is one of the main topics of this work. 

We shall start with justifying our choice of gauge and the quantum picture utilized in the quantization of the system 
atoms+radiation. We then motivate a split of the state space into 0- and 1 -photon states which will be convenient for for- 
mulating the decay of a correlated state in the collectively excited atomic sample. We shall see how the Zw symmetry of the 
configuration can be utilized to diagonalize the effective Hamiltonian of the atomic degrees of freedom on the subspaces car- 
rying the irreducible representations of the symmetry group Z;v; since, for simply-excited states, these subspaces are only one- 
or two-dimensional, the problem of diagonalization then becomes trivial and can be carried out analytically, producing both the 
complex eigenvalues, and the associated eigenvectors, of the effective channel Hamiltonian. We then use this theory to illustrate 
the mechanism behind super- and subradiance; in particular we show that the divergence of level shifts for vanishing interatomic 
distance is due to the Coulomb dipole-dipole interaction only. It will be demonstrated that, for an ensemble with N outer atoms 
in a circular configuration, there are (N — 1) states which are insensitive to the presence of a central atom: The properties of 
the quasi-particle describing the collective excitation do not change when the atom at the center is removed, simply because the 
latter is not occupied. On the other hand, the two states in a circular ensemble which do occupy the central atom are analogous 
to the i-state wave functions in hydrogen-like systems: Just as the s states carry angular momentum quantum number / = and 
therefore transform under the identity representation of SO(3), our collective p = states correspond to the identity representa- 
tion of the symmetry group Z;v; and just as the s states are the only ones which are nonvanishing at the center-of-symmetry of the 
Coulomb potential in a hydrogen-like system, so are our p = states the only quasi-particle states which occupy the central atom 
at the center-of-symmetry of the circle. Quantum beats between two extreme p = states are possible, one, in which only outer 
atoms are occupied, the other, in which only the central atom is excited. At certain radii of the circle, the two p = levels cross, 
and the beat frequency vanishes, making the population transfer between the extreme configurations aperiodic. Moreover, we 
show numerically that, by increasing the number of atoms in the circle for fixed radius, we arrive at a domain where the minimal 
spontaneous decay rate in the sample decreases exponentially with the number of atoms in the circle. The associated subradiant 
states of the atomic sample will therefore qualify as single-photon traps. In this examination, a critical inter-atomic distance 
of one-half of the dominant wavelength emerges. Finally, we show that the same critical distance emerges in a linear-chain 
configuration of atoms: here, it signifies a jump-like, almost discontinuous, behaviour of the minimal spontaneous decay rate, 
from close to zero to finite values. 



II. HAMILTONIAN OF THE SYSTEM AND ELECTRIC-DIPOLE PICTURE 



We first formulate the decay of a simply-excited correlated state in a sample of N neutral identical two-level atoms, arranged in 
an arbitrary planar pattern, such that all atomic dipole moments are aligned perpendicularly to the plane. The atoms are labelled 
by A = 1,2,..., N, each atom being localized around a center-of-mass R4 in space. These locations are regarded as fixed in the 
sense that the associated center-of-mass degrees of freedom do not take part in the dynamics; as a consequence, the quantities R4 
are c-numbers, not operators. The atoms are assumed to be identical, having a spatial extent on the order of magnitude of a Bohr 
radius ao, and consist of point charges which are labelled as ^a- With respect to wavelengths associated with optical transitions it 
is legitimate to ignore the spatial variation of the electromagnetic field over the extension ao of each atom, so that we may replace 
the field degrees of freedom Ax(x^ a ,f) by Ax(R a ,f)- This step constitutes the long-wavelength approximation. For reasons of 
consistency we then must also replace the inter-atomic Coulomb energy by its lowest-order multipole approximation, which is 
the dipole-dipole energy Y.a<b^ab '• After canonical quantization, the Goppert-Mayer transformation yields the electric- dipole 
Hamiltonian 

2 2 

# = EE^+EE ^a P + / d\ £/fcD(*)a,(k)ta,(k)- (la) 

A a zm Aa A a< p J s=1 

-E<U»E ± (Ra)=//o + #/ ■ (lb) 

A 

Here the unperturbed Hamiltonian Hq, given in (la), contains the sum over atomic Hamiltonians Hqa including the Coulomb 
interaction y Aa p between the internal constituents of atoms A = l,...,N, but without the Coulomb dipole-dipole interaction 

V A g between atoms A and B, for A ^= B; and the normally-ordered free field energy. Since the charge ensembles A = 1,2,... 
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are assumed to have mutual distances Rab = |Ra — Rb| which are much larger than the typical extension «o of the atomic or 
molecular wavefunctions, wavefunctions belonging to different ensembles do not overlap. As a consequence, atomic operators 
associated with different ensembles commute, 

[xAai,PB$j] =0 , f or A^B . (2) 

The electric-dipole interaction is given in (lb), where Ej_ denotes the transverse electric field operator. The interatomic Coulomb 
dipole-dipole interaction V^ B P seems to be conspicuously absent in (1); however, it is a feature of the Goppert-Mayer transfor- 
mation (and more generally, of the Power-Zienau-Woolley transformation yielding leading to the multipolar Hamiltonian) to 
transform this interaction into a part of the transverse electric field, so that the Coulomb interaction emerges in fully retarded 
form as a part of the level-shift operator on the radiationless subspace of the system. This will be seen in formulae (26 - 28) 
below, thus clarifying that the w/er-atomic Coulomb interaction V^g is certainly contained in the Hamiltonian (1), albeit in a 
nonobvious way. 



III. THE DECAY OF A COLLECTIVE ATOMIC EXCITATION 



The simply-excited uncorrelated states of the sample are the product states \e\,g2, • • • ,gN,0) = 1 1 , 0) , . . ., \gi,g2, •• • ,eN,0) = 
\N,0), where \gA), |<?a) denote the ground and excited level of the Ath atom, and |0) is the radiative vacuum. The states \e) 
are supposed to be electronic excitations, corresponding to the fact that the exchanged photons will have optical wavelengths. 
In the discussion before eq. (2) we pointed out that the typical interatomic distance in the sample should be comparable to an 
optical wavelength, so that electronic wavefunctions belonging to excitations of different atoms will certainly not overlap; as a 
consequence we can refrain from using antisymmetrized electronic states (Slater determinants) and use simple product states to 
describe the atomic sample. 

The simply-excited uncorrelated states \A,0) are coupled to the continuum of one-photon states \gi,g2,-" ,gN,ks) = \G,ks), 
where (ks) is the wave vector and s denotes the polarization of the photon. As mentioned above, the simply-excited correlated 
states will be superpositions of the form 

N 

\C)= J> B |fi,0) , (3) 

B=\ 

where the complex coefficients ca are to be determined from the condition that the states (3) be energy eigenstates of a suitable 
effective Hamiltonian. These states are formally very similar to Frenkel excitons [13-16], e.g. in molecular crystals. 

We are interested in the spontaneous radiative decay of a simply-excited correlated atomic state which is coupled to a contin- 
uum of one-photon states. The dominant contribution to this decay will come from a quasi-resonant single-photon transition so 
that two-photon- or higher-photon-number processes can be expected to play a negligible role, and hence it will be admissible 
to truncate the possible quantum states of the radiation field to one-photon states |ks), and the vacuum |0). The electric-dipole- 
admissible single-photon transitions of our correlated atomic sample are 

\c)^\G,ks) , (4a) 
\ C )^\AB,ks) . (4b) 

where, in the second line (4b), the sample emits a photon and makes a transition to a doubly-excited state \AB). Such processes 
are inherently non-resonant and will be neglected, a step which is usually called the rotating-wave approximation [17]. We 
therefore take into account only transitions of the sample to the common ground state, accompanied by emission of a single 
photon (4a). Thus, the state space of our joint system atoms+radiation is spanned by the simply-excited radiationless states 
\A,0), and the continuum of one-photon states \G,ks), both of which are eigenstates of the unperturbed Hamiltonian Hq. We 
now split the state space into a radiationless subspace, spanned by the simply-excited states \A,Q), and the subspace of one- 
photon states |G,ks); the associated projectors are 

N 

Q= £ \A,0)(A,0\ (5) 

A=l 

and P = f d\ £ s |G,ks) (G,ks|. By construction, Q and P commute with Hq. 

We now formulate the decay of a given simply-excited correlated state \c) into the continuum of one-photon states: Let 
t/(f,0) be the evolution operator associated with the Hamiltonian (1) in the Schrodinger picture, and let us assume that, at time 
t = 0, the atomic sample is in a correlated state, the radiation field is in the vacuum state |0), and no correlations between atoms 
and radiation are present, so that the state vector of the total system is |v[/ (= o) = \ c). We wish to compute the probability that the 
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radiation at time t > is still trapped in the system; in other words, the probability that the system at time t can be found within 
the radiationless g-space, 



P(t)=£(c\uHt,0)\A 1 0)(A,0\U(t,0)\c) . (6) 



N 

L 

A=l 

The evolution operator U(t,0) can be expressed as a Fourier transform over the Green operator, so that, for t > 0, 

(A,0\U(t,0)\c) = - J £^JdEe-i Et ■ (A,0\QG(E+)Q\B,0) c B , (7) 

B= 1 

where we have used eq. (3). The g-space Green operator QG(E+)Q can be computed by standard methods as 

QG(z)Q^q(z-H -QH i Q-QH i P(z-Hq-PH i P)- 1 PH i q] . (8) 



In the present case, the P-space is spanned by one-photon states only, so that PHjP = 0; and similarly, QHjQ = 0. The energy- 
dependent non-Hermitean Hamiltonian in the g-channel 

QX(z)Q = H Q + QH ! P(z-H y l PH I Q (9) 

has left and right eigenvectors 

(p*,z\X(z)=A p (z)(p*,z\ , (10a) 

X(z)\p,z)=A p (z)\p,z) (10b) 

which satisfy a generalized orthonormality relation [18] 

(P*,z\q,z) =8 pg , (10c) 

and which are g-space-complete in the sense that 

G = I>,z}</>*,z| • (10d) 
p 

The projector property, Q 2 = Q, immediately follows from (10c) and (lOd). 

With the help of relations (10) we can write the non-Hermitean channel Hamiltonian in the form 

#(z)=I>,z>A p (z)</>*,z| , (ID 
p 

which is obviously a generalization of the diagonal form of Hermitean Hamiltonians in the associated basis of energy eigenvec- 
tors. The g-space Green operator can be expressed similarly, 

If this is inserted into (7) for z = E + ;'e we obtain an expression for the transition amplitude which is still exact. However, in the 
neighbourhood of the unperturbed initial energy , 

Ei = (N-l)-E g + E e =E G + E eg , (13a) 
Eq=N ■ Eg , (13b) 
Egg — Eg Eg — Tickgg . ( 1 3c) 

the variation of QG(z)Q near z = E + ;'e due to the quasi-resonant behaviour of the denominator can be expected to be much 
stronger than the variation due to the functional dependence z i— > (z); for this reason, M (z) can be replaced by its value H (Ei + ) 
at the dominant energy of the process under consideration. Then (12) becomes 



-1 N 



\p,E,+) (p*,E, 



QG(z)Q^[z-X(E i+ )) . (14) 
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We need the matrix representation of this operator in the uncorrelated basis \A,0). Let us abbreviate 

\C p ) = \p,E i+ ) , (c* p \ = (p*,E i+ \ , A p = A p (E i+ ) , (15) 
then the matrix of (14) in the basis |A,0) is given by 

* (A,0\Cp}(c*\B,0) 
(A,0\QG(E + )Q\B,0) = £ y ' K '' l ' . (16) 
p=l E+-A p 

When (16) is inserted into (7), the Fourier integral can be performed by means of the method of residues, since the eigenvalues 
A p are now energy-independent. The eigenvalues are generally complex, all imaginary parts being negative, 

amA p <0 for all p =1,..., N , (17) 

since all states in the g-space have a finite lifetime and therefore must eventually decay; this carries over into the condition (17). 
The result for the transition amplitude is then 

(A,0|£/(f,0)|c)=£<A,0|Cp)(c;|fl,0) c B e-i A ?> . (18) 

Bp 

We now go back to equation (9) for the non-Hermitean g-channel Hamiltonian. At the dominant energy Ei+ = Ej + ;'e we can 
write 

X{E i+ ) = QH Q + hA{E i )-i^T{E i ) , (19) 

where A(£' 1 ) contains the level shifts in the g-space, while T(Ei) contains the (spontaneous) decay rates. In terms of the 
uncorrelated basis |A,0) we have 

(A,0\X(E+)\B,0)=E i 8 AB + hA AB (E i )-i^r AB {E i ) , (20a) 

where 

UME) = *(Ai {AMH, ^ iG f lHllm , (20b) 
J s *-[ E — Eg — h(£> 

^F AB {E) = n I cPkY, (A,0\Hi\G,ks)(G,ks\ H i\ B ,0) 8{E-E G -h(o) . (20c) 

The initial energy is £, = Ec + nck ef ,, where k eg is defined in (13c). The decay matrix at the initial energy can be computed by 
integrating out the delta function, 

T AB (E i )=TD l (k eg R AB ) , (21a) 

r^-^r , (21b) 
3 / sinX cosX sinX\ 

where F is the spontaneous emission rate of a single atom in a radiative vacuum. The argument of the function D\ is equal to 2% 
times the distance R AB between atoms A and B measured in units of the wavelength X eg = 2n /k eg of the Bohr transition \e) — > \g). 
A plot of the function X D\ (X) is given in Fig. 1. 

[Remark: In the more general case of arbitrarily aligned atoms there are two correlation functions L>i,L>2 emerging in the 
decay matrix rather than just one. This is indicated by our denotation. For parallely aligned atoms, the contribution of D2 
vanishes.] 

The level-shift matrix A AB and the decay matrix F AB are connected by a dispersion relation which, for the case at hand, reads 

CO 
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The energy Eq = N -E g of the unperturbed joint atomic ground state will be negative in general. Using (22) the level shifts can 
be computed using appropriate complex contour integrals, with the result 



A AB {Ei) 



, X — k eg RAB 
sinX cosX 



X 2 X 3 



nX 2 J 



du ■ 



- u (l+u + u 2 ) 
u 2 +X 2 



(23a) 
(23b) 



The integral on the right-hand side of (23b) can be approximated as follows: We replace the exponential e 11 by 1 for < u < 1; 
and by for u > 1. This leads to the following approximation for the function S(X): 



S(X) 



approx 



3 ( cosX sinX cosX 
2 \~X X 2 



1 

X 3 + %X 2 



l-X 1 
X 



arctan — hi 
X 



1, 1+1 
- In =- 

2 X 2 



2l 



(24) 



A plot of the exact function S(X) as given in (23b), and its approximation (24), is given in Fig. 1. 

The spatial dependence of the correlations between the atoms is now fully contained in the functions D\ (X) ("decay") and 
S(X) ("shift"). Both functions tend to zero as X — > °°; this is clear on physical grounds, since X = 27t ^ A ~ Rfl implies that the 
interatomic distance tends to infinity in this limit, in which case all correlations between the atoms must cease to exist. The 
only exception of course is the self-correlation of atom A, as expressed by the fact that the diagonal elements Taa are equal to 
one at each distance. Thus, the decay matrix Tab tends to T ■ 1 as X — > °°, implying that the spontaneous emission rate for any 
simply-excited collective state \c) tends to the single-atom rate Y, as expected. 

Now we study the opposite limit: The decay matrix T A b has a well-defined limit for vanishing interatomic distance, since 



lim Di (X) 



1 



(25) 



and as a consequence, Tab — > T \/A,B in this limit. The reason for this behaviour is as follows: The limit X — > defines 
the small-sample limit, in which all atomic dipoles are within each other's near zone. It is easy to compute that, in this limit, 
their degree of correlation is proportional to the cosine e A • e#, where e A and are unit vectors in the direction of the atomic 
dipoles. Thus, in the small-sample approximation, parallely aligned dipoles are maximally correlated while dipoles oriented 
perpendicularly are not correlated at all. In our sample, all dipoles are parallel, hence the limit (25). 

On the other hand, the function S(X) diverges in the limit X — > 0. The reason for this divergence is different depending on 
whether off-diagonal or diagonal elements of Aab are considered. Let us first study the off-diagonal case: 

Here, Rab ^ 0, and we can rewrite (23a) in the form 



1 



47te /i R, 



/(Rab) 



(26) 



AB 



where /(Rab) is an analytic function of Rab which tends to \ as Rab tends to zero. This last result should be compared with the 
Coulomb dipole-dipole interaction between the two atoms A and B, 



„dip _ &a • dfi - 3 (d A • cab) (dg • Cab) 

AB ~ ^ ~ 



4ke R ab 



(27a) 



For two classical dipoles which are parallel, <1a • = d 2 , and which are aligned perpendicularly to their axis Cab of connection, 
dA • Cab = tl/j • eAB = 0, (27a) takes the simpler form 



V^(Rab) = 



4ne R 3 A 



(27b) 



AB 



It then follows that the level shifts in (26) contain the dipole-dipole Coulomb interaction between the charged ensembles A and 
B, since (27b) is contained as a factor in the off-diagonal matrix element SAaj(£;) describing the level shift of the system due 
to the coupling to the radiative degrees of freedom as well as the dipole-dipole Coulomb interaction between atoms A and B. So 
we could indeed write 



hA A B(E i )=V^(RAB)-f(RAB) , 

where / now contains the effect of the purely radiative degrees of freedom. As mentioned above, in the limit R A b 
/(Rab) -> \\ eq- (28) then implies that 



Ik Ma(Ei) = \ lim V^(R AB ) 



A^B 



(28) 
we have 

(29) 
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We now turn to discuss the behaviour of the diagonal elements Aaa- These quantities are infinite, since Raa = 0. This 
divergence remains even if only a single atom is considered, and has to do with the inevitable coupling of the atom to the 
longitudinal and transverse degrees of freedom of the radiation field. From (26 - 28) we learn that part of the divergence comes 
from the Coulomb dipole self-energy of atom A. The other contribution to the self-energy is associated with emission and 
reabsorption of virtual photons of the dipole which is an ongoing process even if only one atom is embedded in the vacuum. 
It emerges explicitly as a radiative dipole self-energy (an infinite c-number term), when the Goppert-Mayer transformation is 
performed on the Standard Hamiltonian. These processes are usually summarized by saying that the energy levels, together with 
the states, undergo radiative corrections. 

As a consequence, the unperturbed states are not really the true physical states; rather, they are abstract constructions whose 
coupling to the continuum of radiation modes shifts their unperturbed energy by an infinite, but unobservable amount, an effect 
which is well-known in Quantum Field Theory [19-24]. This amount must be absorbed into the definition of the unperturbed 
energy, so that the energy of the initial state |\|/ f =o) = \c) becomes renormalized 

Ei — ► Ei+ lim = % , (30) 

where Ej is now assumed to have a finite value, namely the value unperturbed by the presence of N — 1 other atoms; as a 
consequence, £, must be assumed to have been infinite in the first place. In turn, we now must subtract the same infinite quantity 
from the diagonal elements of the level-shift matrix, which amounts to a redefinition 

A AB m^~A AB m^{ > A /J B (3D 

As a consequence, the only observable level-shifts are now those due to interatomic interactions. 
Using the functions D\ and S, the matrix element (20a) can be written as 



(A, 0| |£,0) 



fir 

' 2 



S(k eg R AB ) + iD\ (k eg R AB ) 



Ei-i r § , A = B 



E^ab^—Hab , (32) 



where the matrix elements are determined by the function M(X) = S(X) +iD\(X), which can be determined from (21c) and 
(24) to be 

3 iX / I . i 1 \ 3 fl-X 2 1.1, l+X 2 



M(X)=S(X) + iDl (X),-e« ^_ + _-_j + _^_« rcto _ + i + - hl _j . (33) 

Evidently, if we have found the left/right eigenvectors of then we have diagonalized the whole channel Hamiltonian, for yi p 
is an eigenvalue of if and only if Ej — 4f n p is an eigenvalue of Ej 1# — ^ ^ . It follows from eq. (11) that the matrix of the 
channel Hamiltonian M (£,) in the basis |A,0) can be written as 

(A,OK(£0|B,0)=£(A,0|c p )|£ ( --S^/i p |(c;|B,0) . (34) 

We must be aware that £, has been renormalized by an infinite amount and is now to be regarded as the finite, physically 
meaningful energy of the simply-excited atomic sample when all interatomic distances are infinite. A comparison of (34) with 
(20a) now shows that the quantities — ^/j p are the right eigenvalues of the matrix 

M AB (E ( -)-*'2 r Afl(£») . (35) 

It is therefore appropriate to denote the real and imaginary parts of these eigenvalues according to 

nr ( \ ~ h 

- — MEHe/ip + Dm n P = HA P - i-T p , (36) 



and hence 



r 

A p = --fHeyu P , (37a) 
r p = Tlmnp , (37b) 



8 



where A a now denotes the renormalized level-shifts. A comparison of (34) with (1 1), on the other hand, shows that the quantities 
in curly brackets in (34) are equal to A p . From (34, 35, 36) it then follows that A p should be decomposed according to 

A p = E i + HA p -i^r p , r„ >0 , (38) 
where all imaginary parts F p must be positive due to (17). 




FIG. 1: Plots of the correlation functions D\(X) and S(X) which are responsible for spontaneous decay ("D") and level shifts ("S"), respec- 
tively. D\ becomes maximal at zero interatomic distance while S tends to minus infinity there, due to Coulomb- and radiative self-energy of the 
atomic dipole. Both functions tend to zero as X — > °°, expressing the fact that any kind of correlation must cease to exist for infinite distance. 
The inset shows a comparison between the exact shift function S(X), based on eq. (23b), and the approximation S a pprox, based on eq. (24), 
respectively. 



We now insert (38) into equation (18) for the transition amplitude and compute our target quantity, namely the probability 
P(t) for the atomic sample at time t > to be found in the radiationless g-space, as defined in (6); the result is 

P(t) = E c k P B 2 B l c Bl , (39a) 

B t ,B 2 

N _ 

Pb 2Bi = E (B 2 ,0\ c; 2 )(C P2 \A,0) e'( A P2-^)' 

A,pi,P2=l 

xe -l( r >i+ r >2> (A,0\ C PI ) (c* pi \B h 0) . (39b) 

For the given radius r, number of atoms N, and configuration (a) or (b), let A m ; n be the eigenvalue with the smallest decay 
rate, 

r min <r ; , for all p =1,...,JV . (40) 

Then the associated correlated state 

\C)= n r \ r J ^n) (41) 
V \Lmin| Lmin/ 

has the longest lifetime with respect to spontaneous decay [saying nothing about stability against environmental perturbations], 
and hence is a candidate for a single-photon trap. 



IV. CYCLIC SYMMETRY OF THE CIRCULAR ATOMIC CONFIGURATIONS 



In the remaining part of this paper we shall apply the above theory mainly to circular configurations of atoms, containing 
atoms along the perimeter of the circle, and an optional atom at the center. The system with [without] central atom is referred to 
respectively as configuration a) [ b)], see Fig. 2. In the last section we briefly give some results about configuration c), involving 
jV atoms arranged along a straight line. 
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Now let us study the circular configurations: We shall construct the eigenvectors | C p ) of H (£,•) explicitly by group-theoretical 
means, taking advantage of the fact that the system has a cyclic symmetry group 



(42) 



where the generator T is realized by a (passive) rotation at the angle 2k /N about the symmetry axis of the circle, and N is the 
number of outer atoms along the perimeter of the circle. This holds for both configurations a) and b). All cyclic groups of the 
same order N are isomorphic to the group of integers Z# = {0, 1, . . . , Af — 1} with group operation A o B = A + B mod N, unit 
element 0, and inverses A -1 = —A mod N. The set of all unitary irreducible representations of Z# is given by 



/ 2nipA \ 



1,2,..., N 



where p = 0,l, 



p>(7-)=ex P ^-j 
,N—1. The atoms are located at the center and along the perimeter of a circle with radius r such that 

R z = (o,o) , 

2jc(A-1) . 2jc(A-1) 
R A = r [ cos — , sin 



(43) 



N 



N 



(44) 



forA= 1,...,N. 





N+l atoms 


N atoms 
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N atoms 





FIG. 2: The configurations of N(+l) atoms which are examined: In configuration (a), N outer atoms form a regular A'-polygon with radius r, 
plus an atom at the center. Configuration (b) is the same as a), but without the central element. Configuration (c) has N atoms in a linear chain 
with next-neighbour distance d. 

As mentioned above, the basic symmetry operation in our system is the shift T, corresponding to a passive rotation T c about 
the symmetry axis by an angle 2k/ N. T c transforms the coordinates Ra of the outer atoms and the central atom R z according to 



T C R A 

TcRz 



Ra-i 



N+l = l modN 



(45) 



Here R denotes the triple of coordinates with respect to an orthonormal basis which undergoes rotation, and not the invariant 
vector. Next we need to construct the action of the unitary operator T associated with the coordinate transformation T c given in 
(45) on the state vectors of the system. T operates on the atomic degrees of freedom only, so we can ignore the photon states for 
the time being. It is sufficient to specify the action of T on product states; superscripts denote the label of the atoms: 



l<)< 



OK) • 



(46) 



This action derives from eq. (45), and preserves the subspaces of any given degree of excitation. On the single-excitation 
subspace, T acts according to 



r|A) = A-l) , T\z) = \z) , 

and hence has the matrix elements 

T AB = (A\T\B)=8 A , B - 1 . 
The inverse is given by the transposed matrix, T^ 1 = T T , expressing the fact that T is a unitary operator. 



(47) 



(48) 
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Eqs. (47) are consistent with the transformation behaviour of a scalar wave function: Suppose that V|/(x) = (x|\|/) is a wave 
function in the coordinate system x, and T c x = x' is a coordinate transformation. The wave function is a scalar under the 
associated unitary transformation T on the state space if, in the new coordinates x', the new state vector |\|/) with wave function 
V|/(x') = (T c x\Ty) satisfies 



and as a consequence 



V(x')=¥(x) 



T\x) = \T c x) 



(49a) 



(49b) 



This should be compared with the case of the simply-excited correlated states at hand: Here we can consider the sites R z ,Ra 
of the atoms as discrete locations analogous to the parameters x in (49b). The simply-excited state \c) may be regarded as 
describing a scalar quasi-particle (a Frenkel exciton) for which the states \z), \A) are the analogues of the state vectors |x) for an 
ordinary spinless particle. The transformations (45) then imply (47), in full analogy to (49b). 
The channel Hamiltonians of both configurations are invariant under this transformation, 



(50) 



where the matrix ${, was defined in eq. (32). Relations (50) can be proven as follows: The matrix is symmetric by construction, 
and, for configuration (a), has components 



%2z 



%z\ %z2 ■■■ '-KzN \ 



i Hl2 
UlNl %N2 



%1N 



(51) 



such that 



^i = "Kzi = ■ ■ ■ = %zN , 
%a,b = 'Ka+cb+c for \<A,B<N , 

where in (52b) all indices are to be taken modulo N. The relations (52), in turn, follow from the fact that 

R z \=Rz2=-=RzN = r , 
Rab = Ra+c,b+c for 1 <A,B < N , 



(52a) 
(52b) 



(53a) 
(53b) 



and arbitrary C, where r is the radius of the circle. If the matrix representation (48) of the operator T is used on (51) we 
obtain (50). Similarly, for the configuration (b) without central atom, the relevant operators are given by the lower right block 
matrices in eqs. (48) and (51), so that again (50) holds. Thus, we have proven that the generator T of the cyclic group I_n 
commutes with the matrix 'Jt, and hence with the (non-Hermitean) channel Hamiltonian tt(E\), at least on the single-excitation 
subspace; but these arguments can easily be extended to show that T and H commute on all subspaces. Hence we have 



T,SL 







= 



(54) 



Just as in the case of Hermitean operators, the commutativity (54) implies the existence of a common system of eigenvectors 
of 9i and T: Suppose that v is an eigenvector of T, regarded as a complex column vector of dimension N or (N + 1), with 
eigenvalue t , so that Tv = tv. Then 



(55) 



where we have used (54). Eq. (55) says that H preserves all eigenspaces of T, i.e. if T(f) is the subspace corresponding to the 
eigenvalue t, then 



T(f)cT(f) 

On the other hand we must have T A v = t A v, and in particular 



(56) 



tv = T N v = t N v 



(57) 
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from which it follows that 

,=« P (^)=r*(r) (58) 

for some p G {0, l,...,N — 1}, and we see that the eigenvalues of T are just the one-dimensional matrices T p , eq. (43), of 
the irreducible representations of the cyclic group Z#. Accordingly, we can label the eigenspaces of T by the index p of the 
representation as T(t) = T p . 

We see that the eigenvectors v of T in the state space of the atomic system satisfy Tv = T p {T)v and therefore span carrier 
spaces for the irreducible representations (43) of the symmetry group. Since T is unitary, these carrier spaces are orthogonal, and 
their direct sum is the total atomic state space. Since the channel Hamiltonian H preserves the eigenspaces T p of T according 
to (56) we see that a method to simplify the diagonalization of H is given as follows: 

1 . We first determine the eigenspaces T p of T on the state space of the atomic degrees of freedom. These eigenspaces are 
comprised by vectors v each of which transforms as a basis vector for a definite irreducible representation F P (T) of the 
symmetry group Z#, 

Tv = r p (T)v forallveTp . (59) 

2. The channel Hamiltonian M preserves each eigenspace T p , so that we can diagonalize !H on each T p separately. 

Since the eigenspaces T p are smaller in dimension than the original state space we can expect a significant simplification of the 
diagonalization procedure for H . This will be explicitly demonstrated in the following sections. 



V. EIGENSPACES OF THE GENERATOR T OF THE SYMMETRY GROUP 



We now construct the eigenspaces T p of T on the single-excitation subspaces of the atomic systems with and without a central 
atom. These eigenspaces may be constructed by a standard procedure [25] by writing down the projection operators $ p onto T p 
on the state space of atomic degrees of freedom. These projectors have the following property: Given an arbitrary element \c) 
of the (single-excitation) state space, its image <P p \c) transforms like a basis vector for the one-dimensional unitary irreducible 
representation F P (T) defined in (43), 

r A (^|c>)=rp(r A )(^|c>)=exp(^)(^|c>) . (60) 

For the cyclic group 1 N at hand, the projectors are defined by [25] 

where the action of T was defined in (47). In (61a), p ranges between and (N — 1) and covers all irreducible representations. 
The fact that T is unitary implies that the <£ p are Hermitean and have the projector property 

$ p <e 9 = 5 pq <P p ■ (61b) 

The sum over all <P P gives the identity on the atomic state space, 

N-l 

£fp p = l . (61c) 

p=0 

The completeness relation (61c) implies that the state space of the atomic system may be spanned by a basis such that each of 
its members transforms as a basis vector (60) of an irreducible multiplet (In the present case, all multiplets are singlets, since all 
Y p are one-dimensional). 

We now construct such a basis in the single-excitation subspace of the atomic state space which has dimension (N + 1 ) for 
configuration (a) and N for configuration (b). We first treat case (a): We start with an arbitrary correlated state \ c), 

N 

\c)=c z \z,0)+Y,ca\A,0) , (62) 

A=\ 
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and apply the projector £P P ; the result is 



A=Q 



N 



The factor before c z is equal to 8 p o; thus we see that the general form of a normalized basis vector carrying the irreducible 
representation F 1 ' is 



/ 2nipA \ 



fP P| c )=5 p0 c z |z ) 0)+c£ex P ^^j \A 



,0) , N 2 + |c| 2 = i . 



It follows that, for p ^ 0, the eigenspaces of T are one-dimensional, and are spanned by basis vectors 



1 



/ 2nipA \ 



<■■. ) = _£ex P (^pj|A,0) , (C P \C P ) = 1 . 



(64) 



(65) 



This statement is true for both configurations (a) and (b). We emphasize again that, here, N denotes the number of outer atoms. 
c p 

A 



The coefficients c A = (A,0\c p ) of states (65) have the property that 



N- 



= «r 



(66) 



In Fig. 3 we plot the real parts of c p A for some of the ^-states with N = 50 and N = 51 outer atoms. 
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FIG. 3: The real parts $He c p A of the coefficients c p A = (A,0\c p ) for the eigenstates p = 1,2,23,24,25 of a configuration with N = 50 outer 
atoms, and for p = 25 with N = 51 outer atoms. Since p 0, these coefficients are the same for both configurations (a) and (b). 

If p = 0, the associated eigenspace To is one-dimensional for configuration (b), and is spanned by 



1 



(67) 
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Level degeneracy for configurations (a) and (b) 





N=2n+\ 


N=2n 


degenerate 


l,...,n 
2n,. . . ,n+ 1 


p= l,-.,«-l 
2n — 1 , . . . ,n + 1 


non-degenerate 


P=0(±) 


p=0(±),« 



TABLE I: Degeneracy of the eigenvalues for quantum numbers p = 0(±), . . . ,N — 1, for odd and even N, for configurations (a) and (b). 
Vertically stacked p values in the second row are mutually degenerate. 



For configuration (a), the eigenspace To is two-dimensional, with basis vectors | Co) as in (67) and a second basis vector 

\c z ) = \z,0) . (68) 

On account of (61b), basis vectors pertaining to different irreducible representations are automatically orthogonal. For the case 
of the two-dimensional p = subspace of configuration (a) the two basis elements may always be chosen as orthonormal, as we 
have done in (67, 68). 

Thus we have accomplished the decomposition of the space of atomic degrees of freedom into carrier spaces T p for the 
irreducible representations F p of the symmetry group Zn for both configurations (a) and (b). All that remains to be done now is 
to diagonalize the channel Hamiltonian 9i on each of these subspaces separately. 



VI. DIAGONALIZATION OF THE CHANNEL HAMILTONIAN ON CARRIER SPACES OF THE SYMMETRY GROUP 

A. Configuration (b) 

We first discuss configuration (b) without central atom. In this case, each of the eigenspaces of T is one-dimensional, and is 
spanned by states (65, 67). Since H commutes with T, it follows that each of the \C P ), p = 0,...,N—l,is automatically a right 
eigenvector of M , or equivalently, of the matrix . The associated eigenvalue is found to be 

H P = i + j^M(k eg R lA ) cos ( 2%iP ^~ 1} ) , (69) 

which shows that 

p N - p = p p , (70) 
hence some of the eigenvalues are degenerate; the exact result is presented in Table I. On making use of (33) we finally find 



p p =ts(k eg R lA ) cos { ^ %P{A N 1} ) + J>(V?ia) cos ( ^ %P{A N l) 



(71) 



where the real/imaginary parts contains the level shifts and decay rates, respectively. A plot of these quantities for the first four 
eigenvalues for N = 7 outer atoms is given in Fig. 4, where the approximated shift function 5 ap prox based on eq. (24) has been 
used. 



B. Mechanisms for super- and subradiance 



The Figures 3 and 4 give some insight into the mechanism of super- and subradiance. Studying first the limit of vanishing 
radius r — > 0, hence vanishing interatomic distance, we see in Fig. 4 that the p = state is the only one which has a nonvanishing 
decay rate in this limit: All atomic dipoles are aligned parallely in this state [see eq. (67)], and vanishing distance on a length 
scale of wavelength means that all radiation emitted coherently from the sample must interfere completely constructively. Hence, 
the sample radiates faster than a single atom by a factor N, since the decay rate is 



limr„ = o =N-T 



(72) 
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as follows from eqs. (71) and (25). Conversely, the suppression of spontaneous decay for the states with p > is a consequence 
of the fact that the dipoles have alternating orientations, see eq. (65) and Fig. 3, so that, in the small-sample limit r — > 0, roughly 
one-half of the atoms radiate in phase, while the other half has a phase difference of jc; hence 

limr p>0 = . (73) 

r— >0 

This behaviour is exemplified by the plots of ri,r2,r3 in Fig. 4. However, the destructive interference in the (p > 0) states 
is independent of how adjacent dipoles are oriented. For example, in the limit r <C X, the decay rate should not be noticeably 
affected by rearranging the dipoles in different patterns as long as the 50:50 ratio of parallel-antiparallel dipoles is kept fixed. 
What will be affected by such a redistribution is the level shift, since the Coulomb interaction between the dipoles in the sample 
may change towards more attraction or repulsion between the atoms. By this mechanism we can explain the divergent behaviour 
of the level shifts in the small-sample limit: The shift of the p = state always behaves like 

limAp =0 = +°o , (74) 

which arises from the Coulomb repulsion of the parallely aligned dipoles. On the other hand, the level shifts for the (p > 0) states 
depend on the relative number of parallely aligned dipoles in the immediate neighbourhood of a given dipole, or conversely, on 
the degree of balancing the Coulomb repulsion by optimal pairing of antiparallel dipoles. As a consequence, states for which p is 
close to zero always have positive level shifts, since the Coulomb repulsion between adjacent dipoles is badly balanced, as seen 
in the first two plots in Fig. 3, and the behaviour of Ao, Ai in Fig. 4. On the other hand, states for which p is close to N/2 tend to 
have antiparallel orientation between adjacent dipoles, hence the Coulomb interaction is now largely attractive, which explains 
why these states have negative level shifts in the limit r <C X. This is seen in the last four plots in Fig. 3 and the behaviour of 
A2, A3 in Fig. 4. 

The same mechanism clearly also governs the super- or subradiance of the sample with finite interatomic distance. In this 
case the information about the orientation of surrounding dipoles at sites is contained in the transverse electric field, which 
arrives at the site Ri with a retardation IR4 — Ri|/c. Thus, in addition to the phase difference imparted by the coefficients c p A , 
there is another contribution to the phase from the spatial retardation, which accounts for the dependence of the level shifts and 
decay rates on the radius r. Apart from this additional complication, the physical mechanism determining whether a given state 
is super- or subradiant is clearly the same as in the small-sample limit, and can be traced back to the mutual interference of the 
radiation emitted by each atom, arriving at a given site Ri. This radiation is emitted coherently by the atoms on account of the 
fact that the sample occupies a pure collective state. 




0.5 1 1.5 2 2.5 3 

r [X] 



FIG. 4: Relative level shifts A p /T and decay rates T p /T for the states p = 0, 1,2,3 with N = 7 outer atoms and no central atom. From the 
second plot we see that none of the p-states can generally be declared as super- or subradiant; rather, for each state there exist ranges of the 
radius r for which the state \c p ) is maximally super- or maximally subradiant. The physical reason for this dependence lies in the interference 
due to varying phase differences between spatially retarded radiative contributions from all but one atom, at the site of a given atom. 



C. Configuration (a) 

Now we turn to compute the complex energy eigenvalues and eigenvectors for configuration (a) with central atom. For 
p = 1,...,N — 1, the basis vectors carrying irreducible representations of the symmetry group are the same as before, and are 
given in eq. (65). Consequently, they are also right eigenvectors of the matrix 'J{_. The associated eigenvalues /li p turn out to be 
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the same as for configuration (b) and thus are given by formulae (69, 71). This result means that the presence or absence of the 
central atom makes no difference if the system occupies one of the modes \C P ), p = l,...,N— 1, since the central atom is not 
occupied in the states (65). 

On the other hand, the eigenspace To of T corresponding to the p = representation is now two-dimensional and is spanned 
by r° -basis vectors (67) and (68). Since To is preserved by J{ (£,), we now need to diagonalize the channel Hamiltonian in this 
subspace: The representation of the matrix Si in the basis |Co) ,\C Z ) is 









\Cq)-\C z ) 







(75) 



Thus, we must diagonalize a non-Hermitean matrix of the form %J 
Lb, . |Co,±). where temporarily we have defined 

a + i , 1 



a b 
b i 



«fct = 



and 



|Co+) 



1 



VY+c 2 VT+72 

c = 



2 ±- x /(a-i) 2 +4b 2 , 



\Co-) = 
2b 



The diagonalization yields %J |Co,±) = 



(76) 



VT+c 2 



a-i+ ^/(a-i) 2 +4b 2 



(77a) 
(77b) 



These definitions will undergo some refinement, as we shall explain below. The eigenvectors |Co±) corresponding to the column 
vectors (77) then can be written in terms of a complex "mixing angle" angle 8 such that 



\Co+) = cos9 |Co>+ sin9 \c z ) , |Co-) = -sin6 |Co) +cos6 \c z ) 
where tan 8 = c. Such an angle always exists, but is not unique: 

e = ^ini±^ . 

2 ; — c 



(78) 



(79) 



The eigenvectors (77, 78) are not orthogonal, consistent with the fact that the matrix $J is not Hermitean. However, (77) form 
an orthonormal system together with 



(C * + | = cos9 (C |+sin9(<: z | , (c *_| = -sin9 (Co| +cos9 (c z \ . 
The eigenvectors (78) together with the duals (80) now satisfy generalized orthonormality relations according to (10c), 



and completeness in the two-dimensional space To can be expressed as 

|Co+)(c * + | + |Co-)<Co*_| 



I T 



(80) 



(81) 



(82) 



The eigenvalues /liq± can be expressed in terms of the matrix elements of 



1 N 1 



■ A=2 



Y^M{k eg R lA ) 

A=2 



+ 4NM 2 (k eg r) 



(83) 



where the function M = S + iD\ was defined in eq. (33), and r is the radius of the circle. 

The definition of eigenvalues /jo± as given in eqs. (76, 83) is not yet the final one, however. In (76) we have assigned ^0+ to 
the "positive" square root of the complex number (a — i) 2 + 4b 2 , i.e. the square root with positive real part. From (37a) we then 
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FIG. 5: For N = 10 outer atoms, the level shifts Ao± exhibit level crossing twice on every wavelength, for r > X. At the radius r of a crossing, 
the beat frequency (88) vanishes, preventing any oscillatory population transfer between | Co) and \ c z ). The states |Co±) are defined such that, 
as r — > 0, the level shifts Ao± tend to ±°o, respectively. — Plotted are relative quantities Ao±/r and To±/r. For comparison, level shifts and 
decay rates of the state | Co) without central atom are included. 



see that the p$ + so defined always comes with a negative level shift, and therefore the associated real energy £, + Mo+ is always 
smaller than the energy associated with /jo-. This state of affairs would be acceptable as long as the levels would never cross; 
but crossing they do, as can be seen in Fig. 5: For radii greater than a certain lower bound ro, which depends on the number 
N of outer atoms, we see two level-crossings per unit wavelength. For N = 10 this radius is roughly ro ~ X. At each crossing 
we have to reverse the assignment of square roots to eigenvalues in order to obtain smooth eigenvalues. We therefore have to 
redefine po±, and the associated eigenvectors, in order to take account of this reversal at each crossing. The final task is then to 
uniquely determine which of these smooth eigenvalues is to be labelled po+ and /jo-- A look at Fig. 5 shows that, for r < 0.7A,, 
no crossings occur, so that we can uniquely identify the eigenvalues by their behaviour in the limit r — > 0. In this spirit we finally 
define po± to be that eigenvalue whose associated level shift Ao± = — nRe/Jo±/2 tends towards ±oo, respectively. Physically, 
the dipoles in the associated state |Co+) / || Co+| are aligned parallely, which in the limit r — > produces Coulomb repulsion, and 
hence the positive level shift. It follows that the state (p = 0+) is the natural analogue of the (p = 0) state in configuration (b), 
since they both have the same behaviour at small radii. It is also expected to decay much faster than a single atom, an expectation 
which is indeed confirmed by the behaviour of the imaginary part 3m po+ = To+ /T, which tends to^A^+lasr^O. Again this 
follows the pattern of the p = state in configuration (b). A visual comparison between states (p = 0+) and {p = 0) is given in 
Fig. 5 for N = 10 outer atoms. 

On the other hand, in the state | Co-) / || Co- \\, the central atom is now oriented antiparallelly to the common orientation of the 
outer dipoles, and, in the limit r — > 0, is much stronger occupied than the outer atoms, as follows from eq. (78). Hence, as a 
consequence of Coulomb attraction between the outer atoms and the central atom, the energy is shifted towards — °°, and at the 
same time, the system has become extremely stable against spontaneous decay: This is reflected in the fact that, as r — > 0, the 
decay rate To- tends to zero as well. 

As mentioned above, the states with higher quantum numbers p = 1, . . . ,N — 1 are the same as in configuration (b), and 
have the same eigenvalues. The two levels p — (0±) are non-degenerate, except for accidental degeneracy, and also are non- 
degenerate with the (p > 0) levels. As a consequence, the level degeneracy is similar to case (b), and is again expressed in 
Table! 

In a Gedankenexperiment we may think of switching off the coupling of the central atom to the radiation field; in this case 
$Uz = 0, hence 9 = 0, and the energy eigenstates coincide with \ Co) and \c z ). As soon as the central atom "feels" the radiation 
field, the true eigenstates are (non-unitarily) rotated away from this basis. The modulus | tan 9| of the tangens of the mixing angle 
8 may be taken as a measure of the degree of correlation between the two "unperturbed" states | Co) and | C z ), or as a measure of 
the strength of the interaction that couples the two states. Alternatively, we could take the beat frequency co# of the oscillation 
between the two unperturbed states as a correlation measure: 



D. Quantum beats between unperturbed p = states 



As follows from eqs. (78), the true eigenstates |Co±) are in general linear combinations of the "unperturbed" irreducible basis 
vectors |Co), eq. (67), and \c z ), eq. (68); as a consequence, the true eigenstates will give rise to quantum beats between |Co) 
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and \c z ). In order to determine the beat frequency we compute the amplitude for the transition \c z ) — > |Cn) [both states being 
normalized], using eq. (18), 



(Co\U(t,0)\c z ) = (Co\Co+) e-i^' (c * + |C z ) + (C \C -) e'^"-' (c£_\C z ) , 



iA _r I r * 



where 



A ± =Ej- — fio± , 



as follows from eqs. (37, 38). Using (78, 80) and (85) we obtain 

2 



P\c z )^\c Q )(t) = (Co\U(t,0)\c z ) 



sin 8 cos 8 



2 < e - r o+'+ e - r o-'_2e _ 2( r o++ r o-)f cos 



(a 0+ - A - 



(84) 



(85) 



(86) 



The last term shows that the oscillation between the two states occurs at the beat frequency 

a R = A 0+ -Ao- 

The beat frequency can be expressed in terms of the function M, eq. (33), as 



G>fi = ^9te 



Y,M(k eg R u ) 

A=2 



+ 4-NM 2 (k eg r) . 



(87) 



(88) 
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FIG. 6: The beat frequency eq. (88) for three different values = 10,60, 100 of outer atoms, for radii between 2 and 5 wavelengths. At the 
radius of a level crossing Ao+ = Aq_, the beat frequency vanishes. It can be seen that the greater the number of outer atoms N, the greater 
becomes the radius r () beyond which level crossing occurs. 

For a given number of outer atoms, there can exist discrete radii at which the beat frequency G)r vanishes. From formula (87) 
we see that this is the case precisely when the two levels cross, and hence the normalized states \Co+) and |Co,-) are degenerate 
in energy. This is demonstrated in Figs. 5 and 6. 

The beat frequency (Or so computed has to be treated with a grain of salt, however. The reason is that the dynamics in the 
radiationless g _s P ace does not preserve probability flux, since the latter decays into the P-space when occupying the modes of 
outgoing photons. This is reflected in the presence of damped exponentials in formula (86). Depending on the number of atoms 
involved and the radius of the circle, this damping may be so strong, compared to the amplitude of the beat oscillations, that 
the dynamical behaviour effectively becomes aperiodic, i.e., exhibits no discernable oscillations. An example is given by the 
N = 10 plot in Figure 7. 



E. Analogy between p = states and hydrogen-like s states 



It is interesting to note that for the p ^ states, the central atom is unoccupied, irrespective of the radius of the circle, 
or the number of atoms in the configuration. This means that the central atom takes part in the dynamics only in a p = 
state. This is strongly reminiscent of the behaviour of single-particle wavefunctions in a Coulomb potential, such as a spinless 
electron in a hydrogen atom: In this case, the electronic wavefunction vanishes at the origin of the coordinate system, i.e., at 
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FIG. 7: The probability P(t)/\ sin6cos6| 2 for the transition \c z ) — * \Cq), as given in eq. (86), for N = 10,30,50 outer atoms. The quantum 
beats for N = 10 are practically invisible, making the transfer effectively aperiodic. 



the center-of-symmetry of the potential, for all states with orbital angular momentum quantum number greater than zero. On 
the other hand, in the case of our planar atomic system, the circular configurations also have a center-of-symmetry, namely 
the center of the circle; and, as we have remarked earlier, we can interpret the single-excitation correlated states \c p ) / \\C P \\ 
as wavefunctions of a single quasi-particle which is distributed over the set of discrete locations R z ,Ra corresponding to the 
sites of the atoms. Then the amplitudes (A,0\c p ) = c p A play a role analogous to a spatial wavefunction (x|\|/) = V|/(x); and just 
as the hydrogen-like wavefunctions vanish at the origin for angular momentum quantum numbers / ^ [26], so vanish our 
quasi-particle wavefunctions at the central atom for all quantum numbers other than p = 0. In both cases, the associated wave 
functions are isotropic: The s states transform under the identity representation (7 = 0) of SO (3) in the case of hydrogen, and 
states |Co(±)) under the identity representation (p = 0) of 1_ N in the case of our circular configurations. This means that the 
quantum number p is analogous to the angular momentum quantum number I in the central-potential problem; with hindsight, 
this is not surprising, since both quantum numbers p and / are indices which label the unitary irreducible representations of the 
associated symmetry groups Z;v and 5(9(3), and in both cases a rotational symmetry is involved. 



VII. EXPONENTIAL PHOTON TRAPPING IN THE CIRCULAR CONFIGURATION 



In this section we are interested in the photon-trapping capability of maximally subradiant states in the circular configuration, 
for large numbers of atoms in the circle. To this end we choose a fixed radius, increase the number of atoms in the configuration 
gradually, and, for each number TV, compute the decay rate r m ; n of the maximally subradiant state for the given pair (r,N), in 
configuration (b) only. We then expect a more or less monotonic decrease of r m i n as N increases. But what precisely is the 
law governing this decrease? A numerical investigation gives the following result: In Fig. 8(a) we plot the negative logarithm 
— ln(T m j n /r) of the minimal relative decay rate at the radii r = 1, 1.5,2,2.5 A. for increasing numbers of atoms. We see that from 
a certain number TV = N onwards, which depends on the radius, — ln(T m i n /r) increases approximately linearly with N; in the 
figure, we have roughly N = 14 for r = X, N = 20 for r = 1 .5 X, N = 26 for r = 2 X, N = 33 for r = 2.5X. We also see that the 
slope is a function of the radius r. 

We note that the four curves in Fig. 8(a) imply the existence of a common critical interatomic distance: For large N, the 
next-neighbour distance R nn between two atoms on the perimeter of the circle is roughly equal to 2nr/N; if we compute this 
distance for the pairs of values (r,N) as found above we obtain R nn = 0.45,0.47,0.48,0.48 X, respectively. We see that a critical 
distance of R c w 0.5 X presents itself: If, for fixed radius r, the number of atoms in the configuration is increased, the next- 
neighbour distance R nn decreases; as soon as R = R c is reached, the order of magnitude by which spontaneous emission from 
the maximally subradiant state is suppressed becomes approximately proportional to N. Based upon this reasoning we see that, 
for a given radius r, the critical number N of atoms in configuration (b) is given by N — Then, we infer from Fig. 8(a) that, 
approximately, 

r.-r.r'i'K*'-*' , fovN>N , (89) 

where s(r) determines the slopes of the curves in Fig. 8(a); this function decreases monotonically with r. We must have the limit 
s(r) r ~ > °°i 0, because for large radii all correlations between atoms must cease to exist, and hence r m j n — > T in this limit. Let 
X = 1/r denote the lifetime of the excited level in a single two-level atom; then formula (89) tells us that the lifetime X m i n of a 
maximally subradiant state pertaining to r m ; n increases exponentially with the number of atoms, 



^mm = ^-e s{r)(N - fl) , fovN>N 



(90) 
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(a) Circular configuration, without central atom 
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FIG. 8: The negative logarithm of the minimal relative decay rate, being proportional to the order of magnitude of suppression of spontaneous 
decay, as a function of N, for fixed radii of the circle in Fig. 8(a), and for fixed length L of the linear-chain, in Fig. 8(b). In Fig. 8(a) we see 
that for each radius there exists a critical number N of atoms in the configuration beyond which the order of magnitude by which spontaneous 
emission is suppressed is roughly proportional to N. In the Figure, these critical values are roughly at N = 14, 20, 26, 33 atoms. Beyond these 
values, the next-neighbour distance, which is roughly equal to 2%r/N, becomes smaller than the critical distance X/2. - In Fig. 8(b), the same 
analysis is shown for four different lengths of the linear chain. A linear regime is still visible, but both the slopes of the plots, as well as the 
order of magnitude by which spontaneous decay is suppressed, is visibly smaller than in the circular configuration. 
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FIG. 9: The minimal decay rate in the linear chain configuration c), for next-neighbour distances d ranging between and 1.5 A,, and N = 
5, 10, 100 atoms in the chain, respectively. For d < X/2 there is substantial suppression of spontaneous decay even for relatively small numbers 
N of atoms. 



VIII. PHOTON TRAPPING IN THE LINEAR CHAIN CONFIGURATION 



In this last section we study a configuration of N identical two-level systems which are arranged in a linear-chain configuration. 
As before we focus attention on simply-excited states only. This system no longer has a symmetry group, so that we have to resort 
to numerical methods to compute energy eigenvectors and eigenvalues. It turns out that the eigenvectors, i.e., the coefficients 
of the correlated states \c) in the uncorrected basis |A,0), are qualitatively very similar to those seen in Fig. 3 for the circular 
configurations, except for some possible modifications at the boundaries of the configuration. The energy eigenvalues exhibit a 
different behaviour compared to the circular system, though. The most notable difference is the occurrence of a region d e [0,d'], 
d' < X/2, in which spontaneous decay can be strongly suppressed for the state which, for the given distance d, has the minimal 
decay rate r mnl . This is shown in Fig. 9. Furthermore, in the limit N — > °°, the value of d' approaches X/2; and, as soon as 
d = X/2 is approached and exceeded, the minimal decay rate exhibits a jump-like behaviour. This indicates that, even though 
spontaneous decay is suppressed for d < d', this suppression will become increasingly unstable, and susceptible to environmental 
perturbations, as soon as we approach the critical value X/2. In the region below this critical distance, spontaneous decay may 
be suppressed by several orders of magnitude; for example, for the three values of N = 5, 10, 100 as shown in Fig. 9, we have at 
the distance d = 0.25 A. minimal decay rates of r mnl ~ 10~ 2 , 10~ 3 , 10~ 7 , respectively. However, the photon-trapping properties 
of the linear chain are not nearly as pronounced as those of the circular configurations. In Fig. 8(b) we plot the same quantities 
as in Fig. 8(a) but for the linear chain with a prescribed length, fixed at values of L = 0.05,0.1,0.2,0.5 wavelengths. We see that 
a linear regime exists here as well, but the suppression of spontaneous decay is visibly smaller than in the circular configuration, 
as — ln(r m ; n /r) does not exceed a value of 10 in the linear chain, while it reaches close to 40 in the circle. To further exemplify 
this point, let us prescribe a number of, say, N = 40 atoms, and demand a degree of suppression of spontaneous decay of, say, 
— ln(r m i n /r) = 5. We then ask: how small must the next-neighbour distance between atoms in the circle on the one hand, 
and the linear-chain configuration on the other hand, be in order to achieve the prescribed suppression of radiative decay? A 
short computation shows that for the circle we need a next-neighbour distance R nn ~ 2%r/N ~ 0.4A.; while in the linear-chain 
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configuration, d nn ~ 0.005A- is required. In other words, it is much harder to achieve the same degree of photon trapping in the 
linear chain, and indeed, for the values of parameters as given above, the circular configuration performs better than the linear 
chain by a factor 0.4/0.005 = 80 : 1 ! Thus, if our objective is photon trapping, the circular configurations definitely would be 
our first choice. 



We develop the theory of simply-excited correlated states, their level shifts and decay rates, for planar configurations of 
identical two-level systems with parallel dipole moments, and apply the results to circular and linear configurations of atoms. 
For the circular systems, the atomic state space can be decomposed into carrier spaces pertaining to the various irreducible 
representations of the symmetry group Zat of the system. Accordingly, the channel Hamiltonian on the radiationless subspace 
can be diagonalized on each carrier space separately, making an analytic computation of eigenvectors and eigenvalues feasable. 
Each eigenvector can be uniquely labeled by the index p of this representation. For quantum numbers p > the circular 
configuration is insensitive to the presence or absence of a central atom, so that the wavefunction of the associated quasi-particle 
describing the collective excitation of the sample occupies the central atom only in a p = state. It is explained how this 
feature is analogous to the behaviour of hydrogen-like s states in a central potential. The presence of a central atom in the 
circular configurations causes level splitting and -crossing of the p = state, in which case damped quantum beats between two 
"extreme" p = configurations occur. For strong damping, the population transfer between the two extreme configurations is 
effectively aperiodic. Finally, a critical number of atoms, corresponding to a next-neighbour distance of X/2 on the circle, exists, 
beyond which the lifetime of the maximally subradiant state increases exponentially with the number of atoms in the circle. 
The significance of this critical distance is exemplified by the jump-like behaviour of the minimal decay rate in linear-chain 
configurations of atoms with next-neighbour distance close to X/2. It is demonstrated that the photon-trapping capability of the 
circular system is pronouncedly better than that of the linear-chain configurations. 



Here we prove a technical result which is used in the main part of the paper: 

Let Ria denote the distance between atoms 1 and A (i.e., outer atoms only). Let F be any function of this distance, F ~F(R\a). 
Let p be an integer. Then 



IX. SUMMARY 



X. APPENDIX 



A=2 




(91) 



Proof: 

The sum S on the left-hand side of (91) can be written as 




(92) 



The sines are equal to 




(93) 



while the distances satisfy the equations 



Ri,N-A+2=Ru ■ 

If this is inserted into (92) we obtain an expression which is the negative of (91), and as a consequence, S must be zero. 



(94) 
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